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We study numerically a 5D hybrid model which incorporates a split fermion scenario 
and bulk neutrinos. We perform a Monte Carlo analysis of the model in order to find 
the regions in the parameter space allowing for realization of the leptogenesis. We find 
that higher order Yukawa terms must be included in order the model to produce a 
CP violation and net baryon number sufficient for the creation of the observed baryon 
asymmetry of the Universe. 

1 Introduction 

The observed matter-antimatter asymmetry of the Universe suggests that the Standard 
Model (SM) of particle physics is not the final theory. From the premise that the Universe 
was matter-antimatter symmetric above the electroweak phase transition temperature 
0(100 GeV) SM predicts the ratio of present baryon number density over photon num- 
ber density to be ng/n-y ~ 10~ 18 [lj. One the other hand, the Wilkinson Microwave 
Anisotropy Probe (WMAP) [2] measured the result 

— = (6.1 ±0.3) x 10~ 10 , (1.1) 

which means there is a large discrepancy between observations and SM. 

Baryogenesis via leptogenesis is considered to be one of the most appealing scenarios 
of producing matter-antimatter asymmetry in the Universe, [31 B]- Contemporary models 
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extend SM with heavy singlet neutrinos that undergo decay to leptons and antileptons, 
thus yielding net lepton number that is converted to net baryon number via sphaleron 
transitions during electroweak phase transition with the temperature ~ (100 GeV) [5]. 
These heavy neutrinos also account for the observed nonzero masses of SM neutrinos in 
the so called seesaw mechanism [6J. 

There are various ways to extend SM to include heavy neutrinos. One approach is to 
consider SO (10) models where the heavy neutrino mass is close to the GUT scale, while 
the other relies on superstring theories and introduces extra dimensions. In these extra 
dimensional models or brane models, our 4-dimensional Universe, the brane, is immersed 
into larger dimensional space-time, the bulk. The heavy neutrino can propagate in the 
bulk while SM particles cannot, E]. The advantage brought by brane models is the 
possibility to lower the fundamental scale of gravity many orders of magnitude below the 
effective gravity scale, Mp\ and also account for the hierarchies within the SM fermion 
families [9]. 

In a hybrid model introduced in [TU] the bulk neutrino model is incorporated with 
the split fermion scenario [llj. For every brane neutrino one can assign a bulk neutrino, 
which effectively brings flavor dependence to the model. This also results in a great 
number of undetermined parameters, which detracts from the predictive power of the 
model. Making the brane-bulk couplings flavor-neutral will not help either because then 
there is little leeway for the mass matrix to explain observed neutrino oscillations. In 
split fermion scenarios, on the other hand, SM fermions are centered in the brane at 
different locations and their mixings are due to overlapping wave functions. In order 
to reproduce observed neutrino mixings the relative locations of neutrinos in the brane 
are strictly constrained [12]. However, incorporating bulk neutrinos to the split fermion 
models is a way to get round the fine tuning problems. 

Recently we have studied a Hybrid model of neutrinos [10] from the viewpoint of lep- 
togenesis [13]. The effective system consists of a Kaluza-Klein zero mode propagating in 
the bulk and two brane neutrinos. The mass spectrum includes one light mass eigenstate 
that we identify as the electron neutrino of the Standard Model and two heavy states 
that are nearly degenerate. Our conclusion was that lepton and consequently baryon 
number are produced mainly via the mixing loop due to the heavy states. Our result 
for the CP violation parameter is more complicated with respect to earlier studies con- 
cerning the heavy particle mixing contribution to leptogenesis as the light neutrino has 
acquired a mass through the Higgs vacuum expectation value. 
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Because of the complexity of the model in j 13 j . it is not easy to discern patterns on 
how CP violation depends on individual parameters while others are fixed. For this 
reason we elaborate on the allowed six parameters by performing a statistical analysis of 
the parameter space in the present paper. The allowed regions are determined by Monte 
Carlo analysis on the basis that the model consisting of two brane neutrino families 
generates sufficient CP violation for baryon number production and that the lightest 
neutrino mass remains within correct boundaries. We also pinpoint which higher order 
diagrams dominate in the CP violation parameter. 



The model consists of two brane neutrinos z/ 12 and one neutrino \1/ = (ip +1 ip_) T that 
lives in the bulk [TUt [T5] . The brane neutrinos have an expression that is a Gaussian 
with respect to the brane location in the bulk y a 



The extra dimension is assumed to possess orbifold compactification and thus the even 
modes obey ip + (—y) = ip+(y) and the odd modes satisfy ip_(—y) = —ip_(y). 



2 The Model and CP Violation 




(2.1) 




(2.2) 



The five dimensional brane-bulk action reads 
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+vl(x,y)gh(x,y) [^(x, y) + e iS °"ip_(x, y)] }+h.c. 



which after integrating over the extra dimension y becomes 
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After integrating out the heavy Kaluza-Klein modes we are left with the effective 3x3 
mass matrix 



where 
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Erf(^(y a + y^))-Erf(^(y a - y/) )) 



(2.7) 



+e 1 ' 



Erf (^(^ + ^))+ Erf (^(^-^)) 



The physical mass eigenvalues obtained from (12. 6 p comprise one light state mi ~ m a/ j 
and two heavier states with m 2j 3 ~ m. 

Lepton and baryon number are produced in the processes shown in Fig. [T] and the 
corresponding antineutrino production diagrams. Diagrams with mass insertion prior 
to the loop conserve total lepton number and vanish when summed over the final state 
neutrino flavors [T4"t [T5] . Since the decaying heavy states are nearly degenerate, we 
expect the CP asymmetry to arise from the mixing diagrams shown. 
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Figure 1: The relevant Feynman diagrams for the process \2 — * XilH' are shown. The 



tree level diagram due to the decay of \2 to a neutrino and Higgs is in Fig 1 (a) Fig 1 (b) 



and 1(c) depict the mixing diagrams due to the the decay of X2 to a light neutrino and 
Higgs. 



We computed CP violation parameter in [T3] and our result is 
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where 
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The factor 6 is similar in form to a. The contribution of the term proportional to 
sin((6 l 2 — 6%)/2) is negligible because the difference 62 — 63 is suppressed by the large scale 
m and thus the values of the sine function become small. Leaving this part out does not 
affect the allowed parameter regions. 

This result for e deviates from the findings of earlier studies for heavy particle mixing 
and consequent resonant leptogenesis [Hf [T6| FlT\ US] . For one thing, the couplings Ay 
do not satisfy Aj — A^ because of the decompositions of the mass eigenstates. The light 
brane neutrino has acquired a mass and thus mixes with the other brane neutrino and KK 
zero mode in the mass matrix. In contrast, SO (10) motivated models SM have massless 
leptons at the energy scale where the heavy neutrinos decay and thus the neutrino mass 
spectrum consists only of heavy states. 

In the next section, we move forward to perform the numerical analysis which contains 
several physical constraints. First, there is the observed baryon number of the Universe 
which can be derived from the CP violation parameter. Also, the processes producing 
CP violation must have rates smaller than the expansion rate of the Universe. Finally, 
the light neutrino mass must obey certain limits. 

3 Numerical Analysis 

We perform Monte Carlo analyses where we randomize over a set of values for 
(R, yi, 2/2) <^2> M 2 a). The physical constraints consist of the CP violation bound 

— 1.6 x 10~ 6 < £ < —1.6 x 10~ 7 , the condition that the heavy neutrino decay rate 
must be less than the expansion rate of the Universe, r < H and the experimental 
bound on the electron neutrino mass m Ve < 2.2 eV [19]. The washout factor k in Yl = 
Ke/g* is estimated to be 0.01-0.1 [20j-[22j, which gives the above mentioned constraint 

— 1.6 x 10~ 6 < e < —1.6 x 10~ 7 for the CP violation parameter. We have run the 
program so that there are around 20000 points for each parameter that satisfy the above 
conditions so in each of the plots there are around 20000 points. In the rest of this 
section, we present plots of various parameter planes where the allowed regions illustrate 
the essential features of the model. 
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First, we have made a scatter plot in the (R,mi) plane, mi being the light mass, 



shown in Fig. 1(a) The size of the extra dimension R is allowed to vary in the interval 
(10~ 1( , 10 _7 )TeV _1 and applying the conditions mentioned in the previous paragraph 
allow for R values (1(T 15 , Kr 9 )TeV _1 . There are distinct bands that are favored, namely 
R ~ 10 -13 TeV -1 and R ~ 10 -11 TeV -1 . A similar pattern is produced in the (R,yip) 
planes with the same favored values of R. 



In the (R,5i^) planes, see Fig. 2(b) (R,5i) plane, similar patterns are observed as 
in (R,mi) and (R,yip) planes with the distinction that now there is more variation 
in the R values which is reflected in the zigzag-like pattern. Clearly, at 81,2 — nn the 
interval KT^TeV" 1 < R < lO^TeV" 1 is disfavored over small R < KT^TeV" 1 and 
large R > 10 _11 TeV _1 values. At <5 12 — nix/2 the situation reverses and the interval 
lO^TeV" 1 < R < 10- 11 TeV" 1 is favored over extreme R values. 

When we consider the dependence between the light and heavy masses, we yet obtain 
a figure with two distinct bands that are slightly tilted compared to the ones in the 



(R,mi) plane. This graph bears a visual resemblance to Fig. 2(a) since the heavy 
masses 7712,3 ~ R^ 1 ^ 2 . A similar dependence is found between 7713 and mi. 



Fig. 2(a) showcases a distinct pattern where two regions of R values are favored, 
namely R ~ 10 _13 TeV _1 and R ~ 10 _11 TeV _1 . Similar behaviour is presented in Fig. 



2(b) where the two bands form a zigzag-like pattern suggesting that R depends on 



the phase angles Sip. When all parameters save R are fixed, we make the following 
observations as the phase angles 8\p are varied. 

The larger values R ~ 10 _11 TeV _1 are reached when the phase angles 8\p lie in the 
I quadrant and are nearly degenerate, e.g. 81 = tt/12 and 82 = 7r/10. When the values 
are further apart (still both in the I quadrant), the size of the extra dimension R has to 
be of the order (10 -13 — 10~ 12 )TeV~ 1 . If both phase angles 61.2 take values in the second 
and third quadrants, the allowed region for R lies in the vicinity of R ~ 10 _13 TeV _1 . 
As the angles move to IV quadrant, larger R values ~ 10 _11 TeV _1 are again favored. If 
the angles lie in different quadrants, then for instance 5\ G [0,7r/2) and 82 € (7r/2,7r) 
require R ~ 10~ 13 TeV~ 1 . The same applies if 5± G [0,7r/2) and 82 G (tc,3tt/2). As 82 
progresses to the IV quadrant, the larger end of R values are preferred. The larger end 
of R values is also required in the cases where 81 G (vr/2, 7r) and 82 G (37r/2, 2tt) as well 
as with 81 G (7r,37r/2) and 82 G (37r/2,27r). 

Fig. [3] illustrates some of the behavior of e as a function of R when the phase angles 
d>i 2 are varied. The solid curve in Fig. [3] corresponds to the set y~\ = 1.0, 7/2 = 2.0, 



5± = n/12 and 82 = 2tt/3, the dashed curve to the set y\ = 1.0, i/2 = 2.0, 8± = n/12 and 
62 = 47r/3, and the dotted curve to the set y\ = 1.0, 2/2 = 2.0, 81 = n/12 and 62 = tt/2. 

Since two values of R clearly stand out in Fig. [21 we have also studied the behaviour 
of the system when R is fixed to the most favored values R < 10 _13 TeV _1 and R < 



10 -11 TeV _1 . The smaller R value produces Fig. 4(a) which shows how the CP phases 
81 and 62 are restricted. When both phases equal a multiple of tt or tt/2, there is a 
void which indicates these values do not produce the correct magnitude for e. This can 
also be found out when the light neutrino mass m\ is plotted against the phase angle 
5± or 62, resulting in four distinct bands away from 8\^ — kn/2 {k = 1,2,3, ...). When 
R ~ 10 _11 TeV _1 , the voids at half-integer phase angle values (81 ~ mr/2, 82 ~ kn/2) are 



present neither in the (^i,^) plot of Fig. 4(a) nor the (mi, Si^) plots. 

The fixing of R influences the relationship between yi^ and ^i^, as well. In Fig. 



4(b) we see how the phase 8% and the brane location y\ depend on each other at R 



10 _13 TeV _1 . Integer and half-integer multiples of 7r are less favored overall, and in the I 
and IV quadrants negative y\ values are favored and in the II and III quadrants positive 
l/i values are favored. A similar pattern is found for 82, §2)- Plots confronting (82, yi) and 
(81, 2/2) show continuous bands over the y\^ ranges and less points at 8\^ — kn/2, k = 
1, 2, 3, .... If R ~ 10 _11 TeV _1 , the half-integer multiples of tt do not appear as unfavored 
regions in the above mentioned scatter plots. 

The comparison of various terms in (12. 8 p shows that the largest contribution to e 
arises from terms ~ m A A^ in the numerator as the linear terms ~ m 4 Aij undergo 
cancellation and yield a factor that is many orders of magnitude smaller than m^Afj. 
Thus, in contrast to earlier studies of resonant leptogenesis [18J, where the higher order 
Yukawa terms ~ Afj were found to cancel, it is crucial that in our case the higher order 
Yukawa terms are taken into account, despite the fact that one would expect them to 
be unimportant due to their higher order in the Yukawa coupling. Numerically they 
turn out to be of equal magnitude with the terms ~ m^Aij ~ m A (m/v) 2 in our model 
at the allowed R scale. Consequently, the main contribution to e comes from two-loop 
diagrams which are reflected in the numerator terms that go as ~ m^Afj ~ m 4 (m/t>) 4 . 



4 Conclusions 



In summary, we conclude that the size of the extra dimension is restricted to a few 
orders of magnitude away from the Planck scale when the constraints from the observed 
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baryon-antibaryon asymmetry and electron neutrino mass are applied and the higher 
order Yukawa terms are taken into account. The heavy neutrino mass lies in the TeV 
scale. If there were more extra dimensions, the sizes of these dimensions could be larger. 
However, this interesting case would require a separate study as the hierarchy is likely 
to change completely w.r.t. the one extra dimension case |10] . 

As our scenario consists, in addition to two heavy neutrinos, just one light neutrino, 
it is effectively a one flavor model. Therefore it cannot be confronted with phenomena 
involving mixing among the three light neutrino flavors z/ e , u^, z/ T , such as neutrino os- 
cillations and neutrinoless double beta decay. One could extend the model to include 
three light neutrino flavors by adding two extra brane neutrinos. In this case the mass 
spectrum would consist of of three light and two heavier neutrinos. It would be subject 
of a further study to figure out whether the mixing among the light neutrinos can be 
realized consistently with observations in such scenario. It should be noted anyhow, that 
since the heavy neutrinos in such a scenario would be in a super-eV scale and further- 
more unstable, the scenario would not be a viable candidate for a 3+2 model recently 
considered as a possible solution the LNSD and MiniBoone anomalies |23|-[25j. 

Our model differs from earlier studies [T7| [21 [18] in that the light brane neutrino has 
a mass term, implying that the mass eigenstates have also a light neutrino component. 
If only the terms O(Aij) are included, cancellations occur in the CP violation parameter 
e, resulting in a CP violation that is too small to explain the observed baryon-antibaryon 
asymmetry of the Universe. As a consequence, the couplings Aij are different in form 
from those in [T7J [HJ HE] and we do not observe the cancellation of higher order Yukawa 
terms ~ A|- in contrast with the earlier studies. Actually, the situation is quite the 
opposite as in our case the linear terms nearly cancel. Thus, e cannot be considered 
only as an expansion in terms of the Yukawa coupling ~ m 2 /v 2 but also the relative 
magnitudes of m 2 /v 2 x m 2j 3 and \rri2 — m^l have to be considered and higher order 
Yukawa terms need to be included for e to be sufficiently large. 
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Figure 2: Fig. 2(a) shows the scatter plot between R and ni\. Clearly, R ~ 10 11 TeV 1 
and R ~ 10 -13 TeV -1 are favored. Fig. 2(b) depicts the scatter plot between R and 8\. 
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Figure 3: The CP violation parameter e as a function of the size of the extra dimension R, 
with -ei(j/i = 1.0,2/ 2 = 2.0,5! = 7r/12,5 2 = 2tt/3), -£ 2 (j7i = 1.0, j/ a = 2.0,^ = n/12,6 2 = 
An/3), -£ 3 (£i = 1.0,J/ 2 = 2.0,5! = tt/12,5 2 = tt/2). 
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Figure 4: Fig. 4(a) shows the scatter plot between 8\ and S2, where the size of the extra 
dimension is fixed, R ~ 10 _13 TeV _1 . The scatter plot between Si and tji is found in Fig. 



4(b) , also with the same fixed R. 



